Persistent Orbital Degeneracy in Carbon Nanotubes 
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The quantum-mechanical orbitals in carbon nanotubes are doubly degenerate over a large number 
of states in the Coulomb blockade regime. We argue that this experimental observation indicates 
that electrons are reflected without mode mixing at the nanotube-metal contacts. Two electrons 
occupying a pair of degenerate orbitals (a "shelf) are found to form a triplet state starting from 
zero magnetic field. Finally, we observe unexpected low-energy excitations at complete filling of a 
four-electron shell. 
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I. INTRODUCTION 



II. LEVEL DEGENERACY 
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Single- wall carbon nanotubes exhibit a variety of inter- 
esting transport phenomena at low temperatures [1]. At 
weak coupling to the metal electrodes the whole length of 
the nanotube may serve as single Quantum Dot, demon- 
strating a Coulomb blockade behavior at low temper- 
atures [2,3]. As the quality of the nanotube samples 
have improved over time, the single-electron conductance 
peaks in the Coulomb blockade regime have been found 
to group in clusters of four [4-10]. This behavior has 
been attributed to pairs of quantum-mechanical orbitals 
that originate in two subbands of the nanotube electronic 
dispersion and are close in energy, forming four-electron 
"shells". 

In this work, we investigate the shell structure in 
metallic carbon nanotubes. We find that in a large per- 
centage of shells the two orbitals are degenerate with a 
possible level mismatch not exceeding one tenth of the 
shell spacing. The existence of degeneracy in real nan- 
otube devices does not follow trivially from the symmetry 
properties of 2D graphite, as claimed in Refs. [5-8]. In- 
deed, a sizable orbital mismatch was found in most stud- 
ies [4-10]. In this paper, we report on experimental ob- 
servation of the orbital degeneracy over a large number of 
states. Our finding also demonstrates that the electrons 
are reflected without mode mixing at the metal contacts. 
By studying magnetic field dependence, we find that two 
electrons occupying a pair of degenerate orbitals form a 
triplet state starting from zero magnetic field. Finally, we 
observe unexpected low-energy excitations when a four- 
electron shell is completely filled, both in the sequential 
tunneling regime and in the inelastic cotunneling. 

The nanotubes were grown by a CVD method using 
CO as a feedstock gas (the details are described in our 
earlier publication [11]), and Cr/Au contacts were pat- 
terned by e-beam lithography. We present results mea- 
sured on metallic nanotubes. We measure the differential 
conductance at temperatures of 0.3 or 1.2 K by a stan- 
dard AC technique at the excitation level of 10-50 /iV 
RMS. 
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FIG. 1. a) Differential conductance of a 750 nm-long nan- 
otube as a function of V ga te measured at the temperature of 
1.2 K. Inset: schematics of the band structure in a metallic 
nanotube. b) Differential conductance of a 600 nm-long nan- 
otube as a function of V ga te measured at the temperature of 
0.3 K. 



Figure la shows the Coulomb blockaded conductance 
measured in a narrow range of gate voltages on a 750 
nm long nanotube at T = 1.2 K. Most single-electron 
conductance peaks cluster in groups of four, where the 
peaks within each cluster have similar heights [4-10]. 
The spacing between the peaks in each cluster is smaller 
than the spacing between the neighboring clusters. Fig- 
ure lb demonstrates a similar peak clustering measured 
on a 600 nm long nanotube at T = 0.3 K. This behav- 
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ior is further illustrated in Figure 2, where we plot the 
addition energies required to add consecutive electrons 
to the nanotube of Figure la. (The "addition energy" 
[12] is measured as the gate voltage spacing between the 
neighboring conductance peaks multiplied by a factor 
a = Cg a te/Ctotai which we later find to be a = 0.13 
from Figure 4.) 

The exceptional reproducibility of the pattern over a 
wide range of gate voltages and the relatively large nan- 
otube length allows us to trace the positions of hundreds 
of conductance peaks within the same sample. We ex- 
tract the peak positions from the successive gate voltage 
sweeps and keep only the peak positions which coincide 
in different sweeps. By doing this, we get rid of the rel- 
atively rare random offset charge events, which rigidly 
shift the conductance curves in narrow ranges of gate 
voltage. One can see extended ranges of gate voltage 
where three successive addition energies are very close 
and each fourth addition energy is significantly (~ 50%) 
larger than the base line (one of such regions is illustrated 
in Figure 2b). This pattern follows from the peak clus- 
tering observed in Figure 1, where the larger addition 
energies correspond to the spacing between the clusters. 
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FIG. 2. a) Addition energies (spacings between the neigh- 
boring conductance peaks multiplied by a) as a function of 
V g ate, measured on the same sample as in Figure la. b) Same 
in a narrower range of V ga te- Blue squares: addition ener- 
gies corresponding to a complete shell filling (EQ dd )\ green 
circles: addition energy corresponding to a half-filled shell 
(£J ); red diamonds: addition energies of the first and the 
third electrons in a shell (Ef dd , E% dd ). E£ dd points are notice- 
ably higher than the base line, c) Histogram of the addition 
energies. 

Let us for now assume the "Constant Interaction 
Model" [12,13] and ignore the mismatch 8 between the 
orbitals in a shell. Then, a constant charging energy E c 
is required to add an electron to a partially filled shell. 
A larger energy E c + A, where A is the shell spacing, is 
required to add an electron to the next shell. As a result, 
the addition energies within each shell are equal to E c , 
while the addition energy required to put an electron to 



a new shell is E c + A, in agreement with observations 
in Figures 2a and b. Experimentally, we obtain E c w 
3.5 meV and A ss 2 meV. The latter agrees well with 
the expected A = irhvp/L (=2 meV for L = 750 nm). 
Such a simple description may be applicable because the 
exchange interaction J and the excess interaction of two 
electrons occupying the same orbital SU are expected to 
be small compared to the shell spacing A, as discussed in 
the next section (see also Appendix) [14,15]. Both J and 
SU are further reduced by the relatively large diameter 
of our nanotubes (~ 3nm). 

The regular four-electron shell filling is disrupted 
at some gate voltages, where electrons occupy non- 
degenerate orbitals. In such cases, one observes pairs of 
conductance peaks, closely spaced and of similar heights 
[6]. The addition energy of the second electron in a 
pair is smaller than the addition energy of the first elec- 
tron ("even-odd oscillations"), as visible in Figure 2a e.g. 
around V ga te = —4.5 V. It is important, that while the 
orbital degeneracy is absent at some gate voltages, it is 
observed over other extended ranges of gate voltage. 

Apparent clustering of four conductance peaks does 
not yet guarantee the orbital degeneracy. Indeed, the 
clusters are visible even when the two orbitals are not 
degenerate, and the level mismatch S could be as high as 
0.2 to 0.4 of A [4,5,8,9]. Also, a relatively large 5 may 
be masked by the high contact transparency that intro- 
duces lifetime broadening of the peaks [4,7,10]. In our 
measurements, the peaks are narrow and the orbital de- 
generacy is revealed: we find, that for a large fraction of 
the shells 5 is negligible. For example, in Figure 2b, the 
largest level mismatch S is observed at V ga te — —3.34 V 
(the second addition energy in a shell is slightly larger 
than its neighbors). Here 5 = 0.2 meV is equal to 0.1 A. 
The low values of the orbital mismatch, observed over a 
large number of states, point to the existence of a mech- 
anism establishing the orbital degeneracy, rather than an 
accidental alignment of levels. 

The orbital degeneracy survives significant variations 
in A that are likely caused by the disorder potential (no- 
tice variations of the inter-cluster spacings in Figure 2b). 
It is known that long-range potential does not effectively 
introduce back scattering and mode mixing in metallic 
nanotubes [17,18]. At the same time, the orbital degener- 
acy can be easily lifted by structural imperfections in the 
nanotube, such as pentagon- heptagon pairs [15], which 
introduce back scattering and mode mixing [17]. Appar- 
ently, these structural defects are relatively rare in our 
samples. 

Experimental observation of the orbital degeneracy 
makes a non-trivial statement regarding the properties 
of the contacts: evidently, reflections at the contacts do 
not mix the subbands. Let us consider these reflections 
in more details. We surmise that one of the degener- 
ate orbitals originates in the vicinity of point fcrj m the 
nanotube dispersion (inset of Figure 1) and another or- 
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bital originates in the vicinity of point — fco [16]. This 
means than an electron in state fco + 5k and moving to 
the right is reflected from the contact to the state at 
ko — Sk (solid circles in the inset of Figure 1), and not 
to the state —fco — Sk. Similarly for the other orbital, 
electron in —fco + Sk is reflected solely to —fco ~ Sk (open 
circles in the inset of Figure 1). Such processes require an 
electron to change the dispersion branch upon reflection 
(solid and dashed lines in the inset of Figure 1). In fact, 
Ref. [19] suggests that the metal contact perturbs the en- 
ergy spectrum of the nanotube, mixing the two branches 
and opening a gap around ±fc . The resulting spectrum 
at the contacts then look like that of a semiconductor, 
and the electron could adiabatically move from fco + Sk 
to fco — Sk upon reflection. 



An alternative picture (see e.g. a sketch in Figure la 
of [9]), where the states at fco + Sk and fco — Sk belong 
to different orbitals, seems unlikely in our case because 
(i) reflection from fco + Sk to —fco — Sk requires a large 
momentum transfer ~ 2fco, i.e. a sharp nanotube-metal 
interface. If so, the scattering to fco — Sk would be signifi- 
cant, and the two orbitals would mix and split in energy, 
contrary to our observations, (ii) Any possible degen- 
eracy between the orbitals made of the states [fco + Sk, 
— fco — Sk] and [fco — Sk, — ko+Sk] would be only accidental. 
Indeed, in this case the "particle in a box" quantization 
conditions yield for the two types of orbitals the ener- 
gies of hvFi^jr ~ ko) ancl hvp(k Q — ^jt), where N 1<2 are 
integers. The two sets of states coincide in energy only 
if 2fco£/7r equals an integer, a condition that cannot be 
satisfied for a majority of states in a real sample. 



Finally, we analyze the statistics of the addition en- 
ergies (for reviews on this subject see [13,20]). We di- 
vide the addition energies into three categories: the ad- 
dition energy corresponding to a complete shell filling 
(Eq M = E c + A blue squares in Figure 2b), the addition 
energy corresponding to a half-filled shell (E% dd = E c + S, 
green circles), and the addition energies of the first and 
the third electrons in a shell {Ef dd = E§ dd = E c , red di- 
amonds). The statistical distribution of these energies is 
plotted by the corresponding colors in Figure 2c. As ex- 
pected, the overall distribution is clearly bimodal, where 
the main lobe corresponds to the baseline of Figure 2a, 
i.e. to the charging energy. The high energy shoulder of 
the distribution is comprised of E c + 5 and E c + A. The 
distribution of 5 has a significant weight at 5 — and 
for more than 40% of states 8 < 0.1A. We conclude that 
in a significant fraction of the shells the level repulsion 
is suppressed, and the level degeneracy is not accidental. 
Finally, it appears that the distribution of A is broader 
than the distribution for E c (as seen in Figures 2b and 
c). 



III. BEHAVIOR IN MAGNETIC FIELD 
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FIG. 3. a) The nanotube conductance as a function of V ga te 
and B± measured on a 600 nm long nanotube. The con- 
ductance is color coded (bright colors correspond to a higher 
conductance, dark colors correspond to a lower conductance) . 
The inferred spin sequence is indicated by numbers overlaying 
the image. For cluster II, the two middle traces exhibit a kink, 
indicating a singlet-triplet transition, b) Addition energies for 
the four electrons in a shell (similar to shell I in Figure 3a) 
measured as a function of B± on a different 600 nm long nan- 
otube. E± dd and E$ dd stay remarkably constant in magnetic 
field, indicating that two electrons in a row (1st and 2nd; 3rd 
and 4th) are added with the same spin projections. The ad- 
dition energy E% dd corresponding to a half-filled shell grows 
in B± with a slope close to 2/io (dashed line). Inset: conduc- 
tance plot from which the addition energies where extracted. 
The symbols correspond to those in the main panel. 

Figure 3a shows the nanotube conductance as a func- 
tion of the gate voltage and magnetic field perpendic- 
ular to the nanotube axis, Bj_, measured on a 600 nm 
long nanotube. At zero field, we see that the peaks form 
two clusters (I and II) of four electrons, similar to Fig- 
ure 1. The perpendicular field couples primarily to the 
spin of the electrons in the nanotube and not to their 
orbital motion, which allows us to study the spin state 
of the system. In cluster I, we observe that the two bot- 
tom traces move down in magnetic field, and the two 
top traces move up. This behavior indicates that pairs of 
consecutive electrons are added to the nanotube with the 
same spin projection (i.e. spin up or spin down). There- 
fore, two electrons occupying a shell form a triplet state 
[4,10], so that the sequence of the total spin of the nan- 
otube as the electrons fill the shell is — 1/2 — 1 — 1/2 — 0. 
Such a behavior is in contrast with alternating and 1 /2 
states found at zero field in shells with non-degenerate 
levels [5,6,8]. 

In cluster II, the second and the third single-electron 
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traces exhibit kinks around 4 Tesla, a behavior observed 
in [5,8]. Namely, at low B±, the second trace from the 
bottom shifts up and the third trace shifts down with 
the field, indicating that the four electrons in shell II en- 
ter the nanotube with alternating spin projections. In 
higher fields, these directions change: the second peak 
shifts down and the third peak shifts up, just as the traces 
in cluster I. Evidently, in shell II a small level mismatch 
exists at zero field. As a result, at fields below 4 Tesla, 
the two electrons form a singlet state on the lower orbital, 
while at fields above it, the Zeeman energy forces the two 
electrons into a triplet state. From the value of the field 
at the kink location, we find the level mismatch in shell 
II of 8 ~ 0.5 meV (relatively small, compared to the shell 
spacing A w 3 meV in this sample). Any possible level 
mismatch in shell I is at least several times smaller. 

To further quantify the observed behavior, in Figure 
3b we present the addition energies measured as a func- 
tion of B±_ on a different 600 nm nanotube in a shell 
similar to shell I in Figure 3a. The cluster of four peaks 
where the data were measured is shown in the inset. The 
four different symbols mark the valleys from which the 
four addition energies (£$-3) were extracted for the main 
panel of Figure 3b. We see that E^ dd and E% dd (trian- 
gles) stay remarkably constant in magnetic field. At the 
same time, E% (circles) grows steadily with magnetic 
field with a slope close to 2^ (dashed line; g-factor in 
nanotubes is expected to be close to 2). These obser- 
vations prove that the first and the second electron are 
added with the same spin projection (spin- up), while the 
third and the fourth electron are both added spin-down. 
In particular, this supports our statement that the two- 
electron ground state is a triplet down to the lowest mag- 
netic fields. Finally, EQ dd (squares) is significantly larger 
(by approximately A) than Ef dd and steadily decreases 
at a slope close to 2/zo, as expected. 

For degenerate levels (8 = 0), the exchange interaction 
J and excess interaction 8U should show up in E^ dd being 
greater than Ef dd at zero magnetic field (see Appendix) . 
From Figure 3b we have to conclude that the possible J 
and 8U are very small, not exceeding 0.1 A [14,15]. It is 
also noticeable from Figure 3b that Ef dd > E$ dd . (De- 
pending on the particular shell, we find that this inequal- 
ity may be reversed, or the two energies may be equal.) 
This behavior indicates that two electrons occupying the 
first (slightly lower in energy) orbital in a shell have a 
stronger repulsion than the two electrons occupying the 
second orbital (see Appendix). 

IV. LOW-ENERGY EXCITATIONS 

Let us turn to the nanotube conductance at relatively 
large source-drain voltages V s d- Figure 4a shows the con- 
ductance map as a function of the gate voltage V ga te and 
Vsd- The dark regions of suppressed conductance along 



the horizontal axis ( "Coulomb diamonds" ) indicate that 
electron transport is blocked and the number of electrons 
in the nanotube is fixed. The sizes of the Coulomb di- 
amonds follow the same regular pattern as in Figure 1: 
three relatively small diamonds of a similar size follow a 
larger diamond. 
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FIG. 4. a) Differential conductance map as a function of 
Vgate (horizontal axis) and the source-drain bias (vertical 
axis). Each fourth Coulomb diamond is clearly larger than 
its neighbors, b) Zoomed-in view of the 47V and 47V — 1 - 
electron diamonds. A low-lying excitation of the 47V electron 
state is shown by the dashed line. Its energy 8E is clearly 
smaller than the shell spacing A. White horizontal bars indi- 
cate the co-tunneling thresholds. 

Figure 4b shows the zoomed-in view of the 47V — 1 and 
47V electron diamonds in Figure 4a. The lines bordering 
the diamonds correspond to the resonances between the 
energy level in the nanotube and the Fermi energy of ei- 
ther lead. We use them to extract the coefficient a [12]. 
Similarly, the inclined lines running parallel to the dia- 
mond boundaries correspond to the resonances between 
the Fermi energy at the contacts and excited states of the 
nanotube in the sequential tunneling regime. The hori- 
zontal lines visible inside the Coulomb diamonds (marked 
by white dashes in Figure 4b) correspond to the inelas- 
tic cotunneling events, where an excitation is created as 
a result of an electron tunneling through the nanotube 
[21]. 

The states marked by the inclined solid and dashed 
lines in Figure 4b correspond to the ground and the ex- 
cited states of the nanotube with 47V electrons (the top- 
most occupied four-electron shell is completely full). We 
find that these states are separated by an energy 8E of 
only ~ 0.7 meV, several times smaller than the shell split- 
ting A extracted the ground state energies in Figure 2 
(the scale of A can also be gauged from the size differ- 
ence between the AN and 47V — 1 diamonds in Figure 
4). This observation is especially surprising, given that 
the ground state with a completely filled topmost four- 
electron shell should be separated by ~ A from the lowest 
excited state. (This is true up to corrections due to J, 
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SU and 5, which in our case are negligible compared to 
A, see discussion in the preceding sections.) The energy 
of the smallest excitation measured through the inelas- 
tic cotunneling threshold inside a 4iV-electron diamond 
(not shown) coincides with SE extracted in the sequential 
tunneling regime. 

We consider several possible mechanisms to account 
for 5E: quantized ID vibrational modes of the nanotubc, 
features in the density of states of the leads, or collective 
electronic excitations: (i) 8E probably does not originate 
in the quantized phonon modes of the nanotube, as it 
would require a relatively long segment of the nanotube 
(tens of nanometers) to be decoupled from the substrate. 
(ii) The same SE is visible in sequential tunneling in reso- 
nance with both the source and drain electrodes. It seems 
unlikely that both leads would possess a similar feature 
in the density of states. (Hi) Alternatively SE may repre- 
sent an energy of a collective electronic excitation. It was 
recently suggested that in ID GaAs channels at low den- 
sities the interactions may renormalize the energies of the 
spin excitations to lower values [22]. In our experiment, 
the density of carriers near the contacts may possibly 
be reduced by opening of the local semiconducting gap 
as described above in Section II. In this case, the low- 
density ID gas near the contacts is strongly correlated 
and may sustain the low energy excitations SE similar to 
the prediction of Rcf. [22]. At the same time, the single- 
particle energy A would be determined by the properties 
of the nanotube "bulk" , where interactions are relatively 
less important. Further work is required to identify the 
origin of the low-energy excitations SE. 



V. CONCLUSIONS 

We investigate the degeneracy of orbitals in metallic 
carbon nanotubes. We argue that the degeneracy does 
not follow directly from the presence of two subbands in 
the electronic dispersion of a nanotubc. Indeed, a sizable 
orbital mismatch was found within a shell in prior studies 
[4-10] . We present an explanation of the degeneracy and 
discuss the necessary conditions for its observation. Two 
electrons filling a pair of degenerate levels are found to 
form a triplet state down to lowest magnetic fields. Fi- 
nally, when a four-electron shell is completely filled, the 
large and clearly identifiable single-particle gap allows us 
to observe unexpected low-energy excitations, which may 
serve as an indication of collective electronic modes. 
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APPENDIX A 

In our analysis we use an effective Hamiltonian in the 
spirit of Rcf. [14]. The expression invariant under ro- 
tations of the total spin is derived starting from the 
Hartree-Fock Hamiltonian and using operator identities 
of Ref. [23]. The resultant expression depends on Ni (0, 
1 or 2), the occupation number of orbital i, and on Sj 
(0, | or 0, respectively), the total spin of the electrons 
occupying this orbital: 

E([N it Si}) - e * N * + \ E u H N i N i 

+ \Y. u u N i( N i - 1) - E J ^ & • ^ + ^i) ( A1 ) 

i ijij 

Here is the bare energy of an electron occupying orbital 
i; Uij and are respectively the direct and exchange in- 
teractions between two electrons on orbitals i and j. For 
the practically relevant case of two topmost (partially) 
occupied orbitals a and (3 the expression reduces to 

E(N a ,Np,S tot ) = ]T eiNi + U a pN a Np 

i—a,(3 

+ \Y1 u a N i( N i - !) 

i—a.(3 

-J a[j (S 2 tot -Sl-Sl+ l -N a N fi ) (A2) 

where S tot is the total spin of the electron system, Here 
we disregarded interactions with the lower, completely 
filled orbitals. The energies of the relevant states with 
different electronic occupations (presuming that orbital 
a is slightly lower than (3) are: 
£7(l,0,±) = e a (le), 
£7(2, 0, 0) = 2e a + U aa (2e singlet), 
£7(1, 1, 1) = e a + ep + U a p - J a p (2c triplet), 
£7(1, 1,0) = e a + ep + U a p + J a p (2e singlet), 
£7(2, 1, \) = 2e a + e (j + 2U afj + U aa - J afj (3e), 
£7(2,2,0) = 2e a + 2e p + 4U afj + U aa + U pp -2J aP (4c). 
The addition energies measured in experiment are sec- 
ond differences of these E(N a , Np, S to tai)- If the ground 
state of two electrons is a triplet, 
Ef d = e fj -e a + U afj - J a0 , 
E% dd = e a -ej3 + U aa + J af:} , and 
Ef d = ep-e a + U afj + Upp - U aa - J aP . 
If we assume U aa — Upp and define 5 = ep — e a (level 
mismatch) , U = U a p and SU = U aa — U (excess interac- 
tion), the expressions for addition energies become 
E add = E add = u - J + S, and 
E add = u + 8U + J -5. 
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